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ABSTRACT 

The absorption lines in the soft X-ray spectrum of MCG-6-30-15 are studied using the Re- 
flection Grating Spectrometer data from the 2001 XMM-Newton 320 ks observation. A line 
search of the full time-averaged spectrum reveals 5 1 absorption lines and one emission line. 
The equivalent widths of the lines are measured and the majority of the lines identified. We 
find lines produced by a broad range of charge states for several elements, including almost 
all the charge states of oxygen and iron, suggesting a broad range of ionization parameters 
is present in the warm absorber. The equivalent widths of the lines are broadly consistent 
with the best fitting warm absorber models from lTurner et alJ J2003h . The equivalent widths 
of the absorption lines allow confidence limits on the column density of the species to be de- 
termined. For O VII a column density of lOi^'-^^-lQi^^-SS cm^'^ is found. This column density 
of O VII, when combined with the inferred Fe 1 absorption, is sufficient to explain the drop in 
flux at 0.7 keV as being due to absorption from the warm absorber. Fitting 1 K-edge absorp- 



tion to the spectrum reveals a column of 10 



10 



17.67 



cm of O I, suggesting an Fe:0 



ratio of ^1:2, consistent with the neutral iron being in the form of iron oxide dust. Variability 
is seen in a few absorption lines, but the majority of the absorption features, including the 
prominent absorption edges, stay constant throughout the observation despite variability in 
the continuum flux. 
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1 INTRODUCTION 

MCG-6-30-15 is a bright, nearby Seyfert 1 galaxy that shows many 
of the typical X-ray spectral features found in other Seyfert Is. The 
X-ray spectrum consists of a power-law continuum with F 1.9, 
broad iron Ka emission line and a soft excess (e.g. iFabian et alj 
l2002h . The spectrum of MCG-6-30-15 also shows a sharp drop in 
flux at ~0.7 keV. In the pre-Chandra/XMM-Newton era, this drop 
was interpreted as being due to O VII and O VIII K-shell edges from 
ionized gas in the environment of the AGN ("the warm absorber"; 
lOtani et ali.l996:.Revnolds. 1997 : George et al..l998 ). 

This interpretation of the drop in flux at 0.7 keV was chal- 
lenged bv'Branduar di-Ravmont et all i200l h. using a ^--^120 ks Re- 
flection Grating Spectrometer (RGS) XMM-Newton observation. 
They explain the drop as the blue wing of a broad OVIII emis- 
sion line. The broadening mechanism of the line was taken to be 
the same relativistic blurring that produces a broad iron line at 
~6 keV. Part of the reason for adopting such a model was an ap- 
parent ~ 16,000 kms~^ discrepancy in the position of the drop in 
flux and the energy of the O VII edge (0.74 keV). Smaller drops 
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at --^0.39 and ~0.54 keV were similarly explained as being due to 
emission from C VI and N VII respectively. 

lLeeetai]j200lh countered these claims with yet higher spec- 
tral resolution High Energy Transmission Grating Spectrometer 
(HETGS) obs ervations taken w ith Chandra. This discrepancy was 
interpreted bv lLee et alJ i200ll) to be due to the presence of an Fe I 
L-shell absorption edge at the 0.7 keV drop, which when combined 
with O VII resonance absorption l ines and edge co uld adequately 
explain the drop in flux at 0.7 keV. lLee et alJ i200ll) concluded the 
Fe I could be in the form of dust embedded in the warm absorber 
(the Dusty Warm Absorber, DWA), either in the form of iron ox- 
ides or silicates, whic h would explain the optical observations of 
iRevnoIds et alj j^93) where E(V — B) reddening su ggests strong 
dust extinction is present (see lBallantvne et al]i2003)) for a discus- 
sion of the location of the dust). 

[Sakoetal. l2003h . however, fitted the DWA model of 

|Lee et al. 1 2001) and the Relativisti c Emission Line (REL) model 
of Branduardi-R avmont et aljj200ll) to the XMM-Newton RGS ob- 
servation and claimed that the REL model provides a much better 
fit to the data than the DWA model. Models of absorption by vari- 
ous species were also fitted to the absorption lines seen in the data. 
ISako et all i2003h claimed that the column densities of O VII and 
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O VIII, derived from the observed absorption lines, are insufficient 
to explain the 0.7 keV drop by absorption. 

MCG-6-30-15 was observed for a further three XMM-Newton 
revolutions in 2001, p roviding 320 ks of data when the source was 
at a higher flux level CFabian et al. 2002) than the previous XMM- 
Newton observation reported by Branduardi-Rayrnont et al. (200JJ 
and Sakoetal. (200j). An analysis of the RGS spectrum unam- 
biguously showed the prese nce of Fe I through the absorption fine 
structure at the L-shell edge jTumer et aj20 0A. The 0.7 keV drop 
could also be adequately fitted by Fe I and O VII absorption. Self- 
consistent models of warm absorption were made using the pho- 
toionization code CLOUDY and fitted to the data along with the REL 
model. Both models provided similar quality fits to the data and 
both were capable of explaining the gross features of the spectrum. 
The DWA model was, however, more physically self-consistent, 
since the REL model requires the underlying continuum to be very 
flat below ~ 1 keV and steep above this energy. Conventional mod- 
els of Seyfert 1 emission do not contain such a break. Moreover, 
the predicted equivalent widths of the emission lines o f the REL 
model inLB randuard i-Ravmont et alj <200lh . iSako et all l2003t) and 
iTurner et: al. 1,2003.) were much larger than the largest equivalent 
widths predicted by ionized disc models and als o lacked the iron 
L emission predicted by those same models (e.g. lBallantvne et alj 
[2002; Rozanska et al. 2002). 

MCG-6-30-15 also shows interesting spectral variability: The 
flux in the iron Ka line is observed to stay constant while the 
flux in the continuum varies jLee et alj|200 □Jyaughan & Edelson 
2001;''Lee et al."2002; Shi h. Iwasawa. & Fabian 2002; Fa bian et al] 
2002). If a spectrum from a low flux period is subtracted from a 
spectrum from a high flux peri od then the resultin g spectrum is ade- 
quately fitted by a power-law I Fabian et alj2002l) . These properties 
ledp^abian & Vaughan (2003) to propose a two component model 
for the emission in MCG-6-30-15 consisting of a varying power- 
law component and an almost constant reflect ion dominated com- 
ponent. This differencing method was used bv lTurneretai]<2003h 
to remove the effects of absorption from the 320 ks RGS obser- 
vation and determine the underlying continuum. It shows that the 
spectrum does not contain any strong relativistic emission lines and 
confirms that absorption is responsible for the 0.7 keV drop. 

In this paper we use the absorption lines seen in the RGS spec- 
trum from the 320 ks XMM-Newton observation to measure the pa- 
rameters of individual species that make up the warm absorber. We 
do not assume a continuum, either DWA or REL, but instead use 
a cubic spline fit to determine the flux level around the absorption 
lines. This way the absorption line parameters can be determined 
without any prior assumption being made about which model is 
correct. This information is used to confirm that a DWA model ac- 
counts for the observed spectral features, specifically the drop in 
flux at 0.7 keV. The paper is organised as follows: Section |2| de- 
scribes the method used to find, identify and measure the absorp- 
tion lines as well as the results obtained by this method. In section|4| 
the evidence for the presence and composition of dust is examined. 
In section |5| the variability of the absorption lines is measured. In 
section|5|the results are discussed. 



2 ABSORPTION LINE PROPERTIES 
2.1 Observations and data reduction 

We observed MCG-6-30-15 with XMM-Newton on 2001 July 31 
to 2001 August 5 (revolutions 301, 302 and 303) for ~ 320 ks 



l lFabian et alj2002h . The Observation Data Files (DDEs) were anal- 
ysed using the standard reduction chains provided by the 5.4.1 re- 
lease of the XMM-Newton Science Analysis Software (S AS v5 .4. 1 ) 
and the associated calibration files. The last few ks of each orbit 
data were removed since these periods showed strong background 
flaring. The final RGS spectra were rebinned so at least 20 counts 
were present in each spectral bin and were fitted to trial model s in 
the X-ray spectral fitting package XSPEC Vl I. 2 jAmauJigg^) . 

The previous XMM-Newton observation of MCG-6-30-15 
was taken on 2000 July 10-11 an d consists of ~ 120 ks of data 
iBranduardi-Ravmont et alj|200ll) . It was reduced using the same 
prescription as described above. 



2.2 Line finding 

To determine the presence of any spectral absorption lines, the 
spectra were divided up into a number of regions, each region be- 
ing analysed separately. The regions were chosen so they did not 
span any of the large drops in flux and contained approximately the 
same number of spectral bins (the regions are shown as individual 
panels in Fig.0. The exception to this was the 0.534-0.714 keV 
region which spans the Fe I edge at 0.702 keV. This was done to 
allow the Ovil Is^ — ls5p line to be included in the same re- 
gion as the other OviI Is^ — Isnp lines. In this case a model of 
Fe I absorption was included ^Turner et al. 2003). Within each re- 
gion the spectra were fitted with a power-law. Galactic abso rption 
(A^H = 4.06 X 10^" cm~^: lElvis. Wilkes. & Lockman|[l98ft) and 
an additional multiplicative model component to account for any 
unknown variations in the underlying continuum either from com- 
plexity in the underlying continuum or weak photoelectric edges. 
A cubic spline function with a number of knots evenly separated in 
energ y across the region was chosen for this additional component 
lPressir3[l99i). The number of knots varied between 5 and 10 
depending on the relative size of the region to the FWHM of the 
line spread function. A model absorption line (with Gaussian opac- 
ity profile) was then fitted with fixed energy to the data at regular 
intervals across the region. The separation between points where 
the line was fitted was made to be less than the FWHM of the in- 
strumental response to ensure that no significant lines were missed. 
Once the line had been moved across the region it was fitted with 
its energy free to vary to the point where the presence of the line 
improved the of the fit the most. The significance of the im- 
provement in upon adding this line was assessed using the F-test 
and a 95 per cent confidence limit. If the line was found to be sig- 
nificant then the line position was recorded and the region around 
this line was removed (within 1 instrumental FWHM either side of 
the line centre). The process was repeated until no more points pro- 
duced significant improvements in the fit. The process produced a 
list of 5 1 lines whose positions span the whole RGS passband (see 
Table0. Because of the significance level used above, a few of the 
weaker found lines may be spurious detections. 



2.3 Line parameter determination 

Once the positions of the significant lines had been determined, 
the equivalent widths of these lines were determined. To determine 
the equivalent widths, the regions of the spectra were again indi- 
vidually fitted to a power-law. Galactic absorption and cubic spline 
(and Fe I in the case of the 0.534-0.713 keV region). Then absorp- 
tion line models were placed at the energies found in section |2^ 
and fitted to the data. 
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It was decided to use a simple notch model for the line model 
rather than the more complicated Gaussian opacity profile model 
used in Section l2!2l A notch model is an absorption model in which 
the transmission is zero within the line width and unity outside. The 
XSPEC fitting package convolves this absorption model with the in- 
strumental response before fitting it to the data. It was found that 
this notch model could be adequately fitted to the data for all but 
one of the line positions^. This suggests that the widths of the ab- 
sorption lines are less than the width of the instrumental line spread 
function. The notch model has the advantage that its width gives a 
direct measurement of the equivalent width of the absorption line 
and so error determination on the equivalent width is more straight 
forward. It also provides a good approximate model for an absorp- 
tion line lying on the saturated (flat) part of the curve of growth. 
The accuracy of this model at reproducing absorption lines that 
have been convolved with the instrumental line spread function was 
tested. The model was found to accurately reproduce the profile 
shape for a broad range of line parameters. 

The best fits to the regions have reduced values of ~ 1.1. 
The best fitting continuum model does not represent the true emis- 
sion from the central source as the amount of continuum absorption 
from photoelectric edges at any given point in the spectrum is un- 
known. Also the details of reflection spectra are highly uncertain 
and any attempt to introduce a self-consistent model of the contin- 
uum will result in the continuum level being at the incorrect level 
for many of the lines. The purpose of this work is to measure the 
absorption line properties without assuming any physical contin- 
uum, either of the DWA or REL type. The continuum model here 
represents the flux level to either side of the narrow absorption lines 
and allows the equivalent width of those lines to be measured. 

The errors on the equivalent widths of the lines were deter- 
mined using a standard Ax^ = 1-0 criterion. Due to the corru- 
gated nature of the space with respect to line position, a Monte 
Carlo technique was used to determine the line position error and 
avoid problems with the XSPEC error command. The number of 
counts in each spectral bin in the source and background spectra 
were perturbed by a Poisson distribution with mean counts equal to 
the source and background count rate respectively. The best fitting 
line model was then fitted to this new spectra and the resulting best 
fitting parameters recorded. This was repeated 1000 times and the 
resulting distribution of parameters used to determine the error. 

The best fitting line models are shown compared to the com- 
bined and fluxed spectrum from all three orbits in Fig. Q and the 
derived line parameters are given in TableQ Also found during the 
fitting was a narrow emission feature at 0.5565 keV with an equiva- 
lent width of 0.777 ± 0.238 eV and is identified as O VII forbidden 
emission. The line finding method of section is only able to 
find the highest equivalent width lines due to the presence of noise 
in the data. Therefore, when the continuum and line model is fitted 
to the data, the continuum component will be below the true level, 
since it is actually fitting to the average level of the continuum plus 
the undetected, low equivalent width lines. Since the model con- 
tinuum is below the true continuum the equivalent widths of the 
model lines will be underpredicted in general. The use of a many 
knot cubic spline to model the continuum means that only those 
missed lines near a found line affect its equivalent width rather 
than an average of all the missed lines in the region. Unresolved 
blends of many lines will also be missed and filled in by the spline 



^ This line (line 43) was fitted with a model line with a Gaussian opacity 
profile. The width of the line was determined to be 1870 ± 220 km . 



continuum . This is true for the iron Unresolved Transition Array 
(UTA, e.g. lSako et al. 2001: Behar et al. 2001) which has the effect 
of lowering the continuum between 0.7 and ~ 0.8 keV. 



2.4 Line identification 

The lines were identified by comparing their positions to the 
energies quoted in several atomic data sources for the po- 
sitio ns of known absorption transitions. Those sources used 
were IVemer et all jl996l) . the NIST^ database and the line lists 



transitions 


iMcLauahlin & Kirbvl Jl998l). 


Behar & Kahn (2002) and Pradhan et alj 



1 2003) were used. Althoueh iVerner et alj <199fll) provides the most 
accurate data (it only includes transition energies that have been 
experimentally verified) it does not include many important transi- 
tions. One example of this is for the highest oscillator stre ngth line 
of Fe XIX which is not present in the lVemer et alJ il99di list. The 
positions for this line in the NIST, CLOUDY, PHOTOION and XSTAR 
line fists are 0.9170, 0.9211, 0.9170 and 0.9030 keV respectively. 
Such variation in energy and the lack of a complete and accurate 
line list makes the correct identification of the observed lines diffi- 
cult. 

In order to identify the observed lines, the transition with the 
largest oscillator strength from each of the species considered was 
compared to the positions of the observed lines. The species con- 
sidered were C V- VI, N V-VII, O I-VIII, Ne II-IX, Si X-XII, S X-XIV, 
ArlX-XVI and FeXVI-XXIV. Since there is such a large variation 
in the theoretical line positions between the atomic data sets, the 
transition positions from each atomic data set were compared si- 
multaneously. If an observed line was close to the transition from 
the data sets and it was reasonable to expect to see it (the transi- 
tion had a high oscillator strength and the element has high cosmic 
abundance or the transition had been seen in other sources) then the 
identification was made. Once the highest oscillator strength line 
of a particular species had been identified, the search proceeded 
for the lines of the same species in order of decreasing oscillator 
strength. If a transition could not be found then any transitions with 
a lower oscillator strength from the same species were not identi- 
fied. The exception to this rule was C V where the Is^ — Isnp with 
n = 3,4,5 lines could be clearly identified and are expected and 
the Is^ — ls2p line would have its location in a confused part of the 
spectrum at 0.3518 keV. No reliable atomic data exists for the inner 
shell transitions of n itrogen although the y are predicted to occur in 
the DWA models of lTumer et alj i2003h . Some of the unidentified 
transitions around 0.42 keV may be due to these transitions. 

Lines from a broad range of ionization states were identified in 
the data. Oxygen and iron states are the most numerous with detec- 
tions of almost all species of oxygen and iron species from Fe XVII 
to Fe XXIV, as well as Fel from the L-edge and Fe II to Fe XVI 
from the UTA iTurner et alj|2003h . At energies below ~0.5 keV 
absorption lines of carbon and nitrogen are significant. A broad 
range of neon species are also detected from Ne V to Ne X includ- 
ing several inner-shell transitions. A careful examination of the data 
shows that several transitions may have been missed because of 
chip gaps and other problems in the instrumental responses. There 



http://physics.mst.gov/cgi-bin/AtData/main_asd 
^ http://www.nublado.org 

^ http://xmm.astro.columbia.edu/photoion/photoion.html 
^ http://heasarc.gsfc.nasa.gov/docs/software/xstar/xstar.html 
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Figure 1. Best fitting line model compared to fluxed and combined RGS spectrum. The lines below the spectrum indicate the positions of chip-gaps and 
other poor quality regions in the detectors. The line labels above the spectrum indicate the positions of the lines for several species. Also labelled is the O vii 
forbidden emission line at 0.5565 keV. 
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Figure 1. Continued 



appears to be a line directly over a chip gap at 0.427 keV which fits 
the position of N VI. Another line appears to exist at 0.5306 keV, 
whose position matches the position of O III. If this line is due 
to O III then every species of oxygen is evident in the spectrum 
suggesting a broad range of ionization parameters is present in the 
warm absorber. T wo velocity zones are identified, as was found in 
ISako et alj J2003'). One has a velocity of 80 ± 260 km and the 
other -1970 ± 160 km s"^ , where the errors are calculated from 
the variance of the individual velocity values for the lines. Only 
two lines (O VII Is^ — ls2p and O VIII Lya) are identified with the 
second, outflowing, zone. Although not all the observed lines have 
identifications some of these may be random Poisson fluctuations 
in the data rather than true lines. 



3 WARM ABSORBER PROPERTIES 
3.1 Confidence region determination 

In Section l23l the equivalent widths of the lines were measured 
in an independent way with no reference to which species might 
be causing those lines or assumptions about the underlying contin- 
uum. Now that those lines have been identified we can use those 
equivalent widths, combined with a knowledge of the atomic pa- 
rameters of those lines, to calculate confidence regions for the tur- 
bulent velocity of the species and its absorbing column density to- 
wards the central source of X-rays. The profile of an absorption 
line which has been Doppler broadened is given by the convolu- 
tion of a Lorentzian profile (caused by radiation damping) and a 
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Table 1. Table showing lines found in the spectra with measured energies, wavelengths, equivalent widths and identifications where they have been made. For 
lines with identifications the rest frame energy and its reference are shown (for t hese identification in br ackets if more than one is present for a single li ne), 
as well as, the velocity relative to the frame of MCG-6-30-15. Reference key: l: IVemer et al] Il996l) . 2: lBehar & NetzeJ |2003) . 3: |Pradhan et alj l2003h . 4: 
[McLaughlin & Kirbv 1 1998). 



No. 


Energy 


Waveleng 


th 


EW 


Identification 


Rest Energy 


Reference 


Velocity 




keV 


A 




eV 




keV 




km s^-*^ 


1 


0.338090 ± 0.000095 


36.672 ± 0, 


,010 


0.372 ± 0.070 










2 


0.364880 ± 0.000076 


33.9794 ± 0, 


.0071 


0.749 ± 0.071 


C VI Lya 


0.367517 


1 


-157 ± 63 


3 


0.36751 ± 0.00012 


33.736 ± 0, 


■ Oil 


0.498 ± 0.059 


C V Is^ — Isip 


0.370923 


1 


461 ± 96 


4 


0.37532 ± 0.00011 


33.034 ± 0, 


.010 


0.407 ± 0.060 


C V Is^ — ls5p 


0.378533 


1 


243 ± 92 


5 


0.37957 ± 0.00017 


32.664 ± 0, 


.014 


0.399 ± 0.092 


C V Is^ — ls6p 


0.382673 


1 


124 ± 133 


6 


0.381635 ± 0.000088 


32.4876 ± 0, 


.0075 


0.547 ± 0.075 


Sxill2s^ — 2s3p 


0.384546 


1 


-37 ± 70 


7 


0.39213 ± 0.00010 


31.6185 ± 0, 


.0082 


0.433 ± 0.071 


ArXII 


0.395210 


1 


34 ± 79 


8 


0.39654 ± 0.00011 


31.2668 ± 0, 


.0085 


0.469 ± 0.087 


SiXII 


0.399752 


1 


108 ± 82 


9 


0.41028 ± 0.000096 


30.2198 ± 0, 


.0071 


0.461 ± 0.072 










10 


0.41735 ± 0.00046 


29.707 ± 0, 


.033 


0.260 ± 0.075 










11 


0.41995 ± 0.00024 


29.524 ± 0, 


.017 


0.321 ± 0.073 










12 


0.431980 ± 0.000081 


28.7014 ± 0, 


.0054 


0.511 ± 0.064 


C VI Ly/3 


0.435562 


1 


164 ± 58 


13 


0.44843 ± 0.00013 


27.6486 ± 0, 


.0081 


0.452 ± 0.071 










14 


0.455880 ± 0.000074 


27.1967 ± 0, 


.0044 


0.602 ± 0.060 


C VI Ly7 


0.459379 


1 


-22 ± 49 


15 


0.466640 ± 0.000093 


26.5696 ± 0, 


.0053 


0.408 ± 0.058 


C VI Ly(5 


0.470402 


1 


94 ± 60 


16 


0.47260 ± 0.00014 


26.2345 ± 0, 


.0076 


0.339 ± 0.062 


C VI Lye 


0.476391 


1 


81 ± 88 


17 


0.48048 ± 0.00015 


25.8042 ± 0, 


.0082 


0.316 ± 0.078 










18 


0.49072 ± 0.00028 


25.266 ± 0, 


.014 


0.284 ± 0.063 


ArXII 


0.495149 


1 


382 ± 170 


19 


0.496510 ± 0.000064 


24.9711 ± 0, 


.0032 


0.806 ± 0.066 


N VII Lya 


0.500324 


1 


-21 ± 39 


20 


0.52284 ± 0.00012 


23.7136 ± 0, 


.0055 


0.594 ± 0.093 


Ol 


0.526790 


4 


-58 ± 70 


21 


0.52704 ± 0.00010 


23.5246 ± 0, 


.0046 


0.934 ± 0.099 


On 


0.531893 


3 


437 ± 59 


22 


0.54035 ± 0.00013 


22.9452 ± 0, 


.0056 


1.67 ± 0.14 


OlV 


0.544507 


3 


-17 ± 74 


23 


0.54959 ± 0.00013 


22.5594 ± 0, 


.0053 


1.123 ± 0.097 


O V 


0.553996 


3 


80 ± 71 


24 


0.558430 ± 0.00127 


22.202 ± 0, 


.051 


0.652 ± 0.086 


O VI KLL 


0.562287 


3 


-253 ± 689 


25 


0.56255 ± 0.00017 


22.0397 ± 0, 


.0065 


0.578 ± 0.091 


O VI KLL 


0.567433 


3 


279 ± 89 


26 


0.56944 ± 0.00013 


21.7730 ± 0, 


.0051 


1.140 ± 0.091 


O VII Is^ — ls2p 


0.573956 


1 


54 ± 71 


27 


0.57358 ± 0.00013 


21.6159 ± 0, 


.0047 


1.006 ± 0.089 


VII Is^ — ls2p 


0.573956 


1 


—2127 ± 65 


28 


0.61606 ± 0.00016 


20.1253 ± 0, 


.0052 


0.62 ± 0.10 


O V 








29 


0.635260 ± 0.000092 


19.5171 ± 0, 


.0028 


0.892 ± 0.067 


O VI 








30 


0.64840 ± 0.00011 


19.1216 ± 0, 


.0031 


0.910 ± 0.077 


O VIII Lyo 


0.653625 


1 


92 ± 49 


31 


0.65252 ± 0.00012 


19.0008 ± 0, 


.0036 


0.950 ± 0.076 


O VIII Lya 


0.653625 


1 


— 1816 ± 57 


32 


0.66070 ± 0.00021 


18.7656 ± 0, 


.0059 


0.876 ± 0.086 


VII Is^ - ls3p 


0.665558 


1 


-119 ± 95 


33 


0.69205 ± 0.00020 


17.9155 ± 0, 


.0053 


0.942 ± 0.087 


VII Is^ — ls4p 


0.697802 


1 


168 ± 89 


34 


0.70651 ± 0.00017 


17.5488 ± 0, 


.0043 


1.11 ± 0.12 


VII Is^ — ls5p 


0.712724 


1 


313 ± 74 


35 


0.76855 ± 0.00035 


16.1322 ± 0, 


.0073 


1.04 ± 0.19 


OviIlLy/3 


0.774623 


1 


46 ± 138 


36 


0.81975 ± 0.00042 


15.1245 ± 0, 


.0077 


0.71 ± 0.22 


FeXVII 








37 


0.82959 ± 0.00034 


14.9452 ± 0, 


.0061 


1.19 ± 0.21 










38 


0.84465 ± 0.00043 


14.6788 ± 0, 


.0075 


1.418 ± 0.021 










39 


0.86513 ± 0.00037 


14.3313 ± 0, 


.0061 


2.12 ± 0.21 


Fe XVIII, (Ne V) 


0.870740 


2 


—379 ± 129 


40 


0.87559 ± 0.00053 


14.1601 ± 0, 


.0086 


1.31 ± 0.20 


Ne VI 


0.884340 


2 


673 ± 183 


41 


0.89000 ± 0.00050 


13.9308 ± 0, 


.0077 


1.54 ± 0.20 


FeXIX, (Ne VII) 


0.897530 


2 


213 ± 167 


42 


0.90043 ± 0.00053 


13.7694 ± 0, 


.0081 


0.76 ±0.22 


NeVIII 


0.908580 


2 


390 ± 178 


43 


0.91337 ± 0.00043 


13.5744 ± 0, 


.0065 


7.63 ± 0.84 


Fe XIX, (Ne IX) 


0.922023 


1 


518 ± 144 


44 


0.95794 ± 0.00050 


12.9427 ±0, 


.0068 


3.46 ± 0.32 


FeXX 








45 


0.97841 ± 0.00063 


12.6720 ±0, 


.0082 


2.06 ± 0.33 










46 


1.00300 ±0.00078 


12.3613 ±0, 


.0096 


2.47 ±0.38 


Fe XXI 








47 


1.0140 ±0.0011 


12.227 ±0, 


.013 


2.88 ± 0.40 


Ne X Lya 


1.021810 


1 


-14 ±321 


48 


1.0486 ± 0.0012 


11.824 ±0, 


.013 


1.49 ± 0.44 


FeXXII 








49 


1.0686 ± 0.0017 


11.602 ±0, 


.018 


1.81 ± 0.45 


Fe XXIII, (Ne IX) 


1.073783 


1 


-869 ± 477 


50 


1.1640 ±0.0032 


10.652 ±0, 


.029 


3.21 ± 0.58 


FeXXIV 








51 


1.34100 ±0.00095 


9.2457 ± 0.0065 


2.31 ± 0.51 


MgXI Is^ - ls2p 


1.352253 


1 


192 ± 214 



Gaussian profile (from Doppler broadening) giving the Voigt pro- 
file. The Doppler broadening is caused by a combination of ther- 
mal broadening and turbulent broadening, with turbulent broad- 
ening being the dominant me chanism for the gas within the sin- 
gle zone warm absorber (INicastro et alJ ll999). The nature of the 
Voigt profile means that for a single line turbulent velocity and col- 



umn density are highly degenerate. To break this degeneracy sev- 
eral other lines from the same species are required or upper limits 
for the equivalent widths if no other lines are positively identified. 
Having more than one line for a species, therefore, allows tighter 
constraints to be placed on the column density and turbulent ve- 
locity than would otherwise be the case. For instance, if the O VII 
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Is^ — ls2p line is considered on its own tlie column density is only 
constrained to lie in the range ~ io^''"''-10^* ** cm~^ and the tur- 
bulent velocity is unconstrained in the 10-350 kms~^ range. The 
expected equivalent width of a line given a particular turbulent ve- 
locity and column density can be found by numerically integrating 
the absorbed area under the normalised Voigt profile. Confidence 
regions for the turbulent velocity and column density of a particu- 
lar species were determined by calculating the expected equivalent 
widths for all the observed lines of that species across a grid of 
turbulent velocity and column density values. These values were 
then compared to the observed equivalent widths and errors cal- 
culated in section I23I to produce a grid of values. By drawing 
contours at particular values of confidence regions could be cal- 
culat ed. The atomic data used to calculate the Voigt profiles was 
from 'Verner et alj il996t) . the same source of atomic data used by 
ISakoetaLaOOl . 

This method can be checked by replacing the notch models 
with Voigt profiles for the lines of interest (constraining the column 
density and turbulent velocity to be the same for all of the lines) and 
fitting to the data. Confidence regions are then obtained by calcu- 
lating the Ax^ from the best fit for a grid of column density and 
turbulent velocity. 



3.2 Results 

3.2.1 Oxygen 

O VII is a prime candidate for analysis since four lines in the first 
velocity zone are identified (Is^ — Isnp, n — 2,3, 4, 5) in a region 
of the spectrum with high signal to noise. The confidence region 
for turbulent velocity and column density considering all four lines 
is shown in Fig.|2| The column density of O VII in the first veloc- 
ity zone is shown to lie in the region ig^^ '^^-lQ^*'^* cm~'^ and the 
turbulent velocity in the range 86-99 km with la. The implied 
optical depth of the O VII edge is 0.57 (taking the lower la bound 
of the column density) suggesting that O VII absorption is responsi- 
ble for the majority of the drop at 0.7 keV as predicted by the DWA 
model. This column density places the Is^ — Isnp, n = 2, 3, 4, 5 
lines on the saturated part of the curve of growth. The result was 
confirmed by fitting Voigt profiles d irectly to t he dat a (see section 
13. II and Fig. |3}. As mentioned by Sako et al. (2003) the presence 
of Ovil forbidden emission suggests that the Ovil Is^ — ls2p 
may be partially filled in b y resonance emiss ion. To take this into 
account the prescription of ISako et alj i2003h was used where the 
ratio of forbidden to resonance emission is fixed at a value of 3 
for a recombination-dominated plasma. One other issue is that the 
Is^ — ls5p transition lies over the position of the Pel L3 edge 
for Galactic absorption. By calculating the equivalent width of the 
line-like feature at the edge energy of Fe I this can be taken into 
account. By taking the limit of the confidence regions from all the 
various combinations of corrections a more robust confidence re- 
gion is obtained. This produces a Ict confidence limit on the col- 
umn density of lois ^S-lQi^ se cm"^ and turbulent velocity of 79- 
lOOkms"^ . If a similar analysis is performed for O VIII then la 
confidence limits for column density of lO^^'^^-lO^*'^^ cm~^ and 
for turbulent velocity of 46-1 1 1 km s^^ are obtained. 

The O VII colum n density lim it obtained above is much larger 
than that obtained by ISako et alJ j2003) using the earlier XMM- 
Newton observation. To investigate this issue the 2000 XMM- 
Newton observation was analysed using the method described 
above. The best fitting model from the 2001 data was used but the 
line and continuum parameters were allowed to vary. Equivalent 



05 




I , , . . I , , , . I . , , , I . , , 
^ 50 100 150 200 

Line Width (b parameter; km s"^) 

Figure 2. Confidence region for column density and turbulent velocity for 
OviI using measured equivalent widths of the OviI Is^ — Isnp, n = 
2, 3, 4, 5 absorption lines. 




50 100 150 200 

Line Width (b parameter; km s"^) 



Figure 3. Confidence region for column density and turbulent velocity for 
OviI determined by fitting Voigt profiles to the OviI Is^ — Isnp, n = 
2, 3, 4, 5 absorption lines. 



widths and errors for the O VII Is^ — Isnp, n = 2, 3, 4, 5 lines 
were measured as above and used to determine a confidence region 
for column density and turbulent velocity for the 2000 data. This is 
shown in Fig.|4]and is compared to the limits found bv lSako et alj 
(2003) and in this paper for the 2001 data. The contours overlap 
the 2001 contours (which are much smaller because of the higher 
signal-to-noise of the 2001 observation) but do not overlap the con- 
tours found by S ako et alJi2003lh . An examination of the equivalent 
width obtained for the O VII Is'^ — ls2p line reveals a significant 
di fference between th at found here (1.14 eV, 44 m A) and that found 
bv lSako et ail <2003h (0.48 eV, 18 mA). We are unable to account 
for this difference. The region around the O VII Is^ — ls2p line 
is shown in Fig.|5| Also shown is the best fitting model from Sec- 
tion l2.3l as well as this model with the O VII l.s^ — ls2p line with 
the equivalent width expected usin g the column densi ty and turbu- 
lent velocity of O VII determined in ISakoetalJ j2003'i. This equiv- 
alent width also includes the expected response emission infilling. 
Clearly the lSako et alj <2003h values of column density and turbu- 
lent velocity under-predict the equivalent width of the line. 
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50 100 150 aOO 350 

Line ffidth (b parameter: km s"^) 

Figure 4. Confidence regions for column density and turbulent velocity for 
O VII. (black): This paper's analysis using the 2000 data, (solid grey) : This 
paper's analysis using the 2001 data, (dashed ergy^i lSako etaP 120031) anal- 
ysis using the 2000 data. 




o 0.555 0.56 0.565 0.57 0.575 0.58 0.585 0.59 



Energy (keV) 

Figure 5. Detail of spectrum around the O VII Is — 2p line. Overlaid is 
the best fitting model from Section |231 (solid) as well as this best fitting 
model with the equivalent width of O VII Is — 2p for the column density 
and turbulent velocity from lSako et alJ (dashed). Also labelled are 

(7) the O VII forbidden emission Une, (2,3) the O VI KLL lines, (4) the O VII 
Is — 2p line for the first velocity zone and (5) the O VII Is — 2p line for 
the second velocity zone. 



Table 2. Comparison of column densities (A'^x ) and turbulent velo cities 
(Vturb) between this work and DWA models in iTumer et alj i20o3) . Nx 
is given i n units of log cm ~ ^ and Vturb units of kms~^ . The DWA 
models in iTumer et air<2003l) assumed a V^urb of 100 km s . 



Species 




DWA models 




This work 




1 


2 


3 








logTVx 


logAfx 


logATx 


logTYx 


Vturb 




cm~^ 


_2 

cm 


cm~^ 


cm~^ 


kms"-*- 


VII 


18.12 


18.18 


18.19 


18.36-18.86 


79-100 


VIII 


18.12 


18.19 


17.93 


17.42-18.47 


46-111 


Cv 


17.74 


17.82 


17.80 


17.27-19.65 


122-153 


CVl 


17.52 


17.57 


17.57 


17.63-17.84 


48-130 



available for the strongest oscillator strength transition (Is^ — ls2p) 
in Cv, the column density and turbulent velocity could not be 
tightly constrained, la confidence limits were obtained of 10^^'^^- 
lO^^'*'^ cm^^ for column density and 48-130 kms^^ for turbulent 
velocity. 



3.2.3 Nitrogen 

The absorption lines of N VII prove more difficult to understand. 
Only one line is detected (Lya) but by using an upper limit on 
the equivalent width of the undetected Ly/J line a confidence re- 
gion for column density and turbulent velocity can be determined. 
This region, however, does not agree with the region obtained by 
directly fitting a Voigt profile to the Lya line. A possible explana- 
tion for this is that the Lya line is not a single line but made up 
of several lines with a velocity separation greater than the turbulent 
width of the lines, but less than the instrumental resolution. Look- 
ing at the abundance of individual specie s in each zone of the DWA 
models presented in iTumer et alJ i2003h it is apparent that while 
other species have their column densities preferentially in one zone, 
N VII has its column density more evenly distributed between the 
zones. Considering the N VII Lya line to be made up of multiple 
lines from different zones brings the limit on column density into 
reasonable agreement with the predictions of the DWA models. 



3.2.2 Carbon 

Cvi has five identified lines (Is — np,n = 2,3,4,5,6). Calcu- 
lating the confidence region for this species produces a poor fit 
and two separate minima because the n = 3 line has a lower 
equivalent width than the n — 4 line when it should have a larger 
equivalent width due to its higher oscillator strength. There are two 
possibilities to explain this: either the n = 3 is underpredicted (be- 
cause unidentified low equivalent width lines nearby mean the fitted 
continuum is lower than the true continuum) or the n = 4 transition 
is blended with another transition (so increasing its apparent equiv- 
alent width). Since no other transition can be found that fits the 
position of the n — 4 transition the former explanation is favoured. 
Removing the n — 3 line from the analysis results in a single 
minimum in the confidence region with column density lying in 
the range IQ^^ '^^-IQ^^-^^ cm~^ and turbulent velocity in the range 
122-153 km (Icr confidence limit). Because no measurement is 



3.2.4 General 

Table |2| compares the column densities and turbule nt velocities 
derived above with those of the DWA models from iTumer et alJ 
(2003). As can be seen reasonable agreement is found between 
the values empirically derived here and those obtained by fitting 
DWA models to the data. Finally we compare the observed equiva- 
lent widths to the equivalent widths predicted by the DWA models. 
As can be seen from Fig.|6|a reasonable match is found although 
the observed absorption lines are found to be larger than predicted 
by a factor of ~ 1.25. This may be due to a combination of dif- 
ferent species having different abundances than the assumed Solar 
abundance or the gas having a different turbulent velocity than the 
assumed 100 km s^^ . For example Ovil requires an increase in 
abundance of ~ 3 — 5 or an increase in turbulent velocity from 
100to~130 km to account for the ~ 1.25 difference between 
observed and predicted equivalent widths. 
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Figure 6. Compaiison of equivalent widths for absoiption lines between 
those observed and those predicted by the CLOUDY DWA models from 
iTumer et al] <2003ft. Circles indic ate model 1, triangles model 2 and stars 
model 3 from^^^^^^^j200^. Only identified lines with accurate line 
details from R^me^^^IjOQ^T are included. 



4 ABSORPTION BY DUST 

The presence of dust in MCG-6-30-1 5 was first suggested by 
[Reynolds et al. 1 199^ t o explain the larp e Balmer decrement found 
in the optical spectrum. lLee et a L'l2001') found that the drop in flux 
at ~ 0.7 keV had the exact shape of absorption by Fel and that 
the Nh implied by the discontinuity at the Fe I L-edge was consis- 
tent with that determined from E{V — B) reddening. This shape 
has also been seen i n the Galac t ic blac k hole candidate Cygnus X- 1 
iSchulz et alJl200A . lLee et alJ iioOll) attributed this absorption to 
dust consisting of iron oxides or silicates (such as olivine) embed- 
ded in the warm absorber. However the sensitivity of the HETGS 
at the neutral oxygen edge was insufficient to measure the column 
of neutral oxygen accurately. Using RGS data Sako et al. 1 2003) 
produced an upper limit at 90 per cent confidence on Ol column 
density of 10^^ cm~^ using the apparent absence of a drop across 
the O I K-edge at 0.538 keV. Such a column density corresponds to 
a 0.5 per cent drop in flux across the O I edge. 

Using our higher quality RGS data we have measured the drop 
in flux across the neutral O I edge. The 0.5-0.555 keV region was 
fitted to a power-law. Galactic absorption, the line profiles found in 
section |3| and an absorption edge at the Ol edge energy. This pro- 
duces an Ol edge optical depth of 0.16-0.23 (Itj) corresponding 
to a column density of 10^^'^^ - 10^^ "'^ c m~^ . Given t he col - 
umn density of Fe I of lO^'^ '^ cm"'^ found in lTumer et alj i2003h . 
this suggests a neutral Fe:0 ratio of ~1:2. This is consistent with 
the dust being i n the form of iron oxides or silicates as found by 
iLeeetalJiioOll) . 

We repeated the above analysis on the 2000 data. Using 90 per 
cent confidence limits we find an optical depth at the edge of 0. 13- 
0.34, consistent with the result found with the longer observation. 
As well as finding a different optical depth across the edge, the 90 
per cent confidence limit range found here is an order of magnitude 
larger than that found bv lSako et alJ j20o3) using the same data. 



5 LINE VARIABILITY 

5.1 Short timescale variations 

To test for short term variability in the warm absorber the observa- 
tion was split into two sets; one when the flux level was below a 
certain level and one when it was above that level. The level was 
set so an equal number of counts was present in each set. The level 
was found to be 16.6 cts for the 0.2-2 keV EPIC pn light-curve. 
The data from each set were then reduced using the method de- 
scribed in section ITTI The best fitting continuum and line model 
from section |3| was fitted to each set of spectra separately and the 
errors for the best fitting parameters calculated as before. No sig- 
nificant change in position of the lines was found. Only 4 lines 
have a ratio of the equivalent widths from the high flux period to 
the low flux period which individually differs significantly from 1 
(greater than 2a; lines 20, 21, 36 and 46 corresponding to O I, O II, 
Fe XVII and Fe XXI with significances of 3. Oct, 3.0a, 4.Sa and 2.0a 
respectively). Given the large number of lines found, it is likely 
that random fluctuations would give false variability detections for 
one or two lines at a 2a threshold (the false alarm probability for 
51 lines at 2a individually is ~0.4). Since three of the lines have 
higher detection significances than this, the variability in these lines 
is likely to be real. Individually all other lines have a ratio that is 
consistent with 1. The variability in line strength can be used to 
estimate an upper limit for the distance between the warm absorb- 
ing region and the central source. Assuming that the increase in 
equivalent width is produced by recombination of more highly ion- 
ized species (brought about by the drop in continuum flux) then 
the recombination timescale. Tree ^ Tots, the timescale over which 
the spectrum is observed to change. Tots is taken to be the average 
length of bright and dim time intervals, which is ~ 2500 s. From 
■Shull & van Steenberg 1 1982) 

^rad-^ ... 
Tree = (1) 

n 

where T is the temperature of the gas, n is the density and A^ad 
and Xrad species dependant parameters. Using t he definition of the 
ionization parameter used in lTumeretalJj2003l) 



then 



AradS.R'T^ 



Consequently, if Tree ^ fobs then 



R<. 



TobsL; 



2-10 Arad 



(2) 



(3) 



(4) 



The Ol and Oil lines are produced solely in the low ioniza- 
tion zone. T aking the appropriate values for ^, L2-10, T from 
iTumer et al .' 1 2003 > gives an upper limit of the radius of this zone 
as 10^^'^ cm. FeXVII and FeXXI are found predominately in the 
highly ionized zones and using the appropriate values for this zone 
constrains it to lie within 10^'^'" cm of the central source. 



5.2 Long timescale variations 

To test the longer timescale variability in the warm absorber, the 
data from a previous observation of MCG-6-30-1 5 with XMM- 
Newton was extracted using the same process as described in Sec- 
tion 12.11 The line and continuum model of Section |3] was fit- 
ted to the data and the resulting parameters compared to the val- 
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Figure 7. Detail of RGS spectra for high flux period {upper) and low flux 
period (lower). The dashed lines show the position of the O I and O II ab- 
sorption lines. Their equivalent width is clearly higher during the low flux 
period. 

ues obtained from the 2001 data analysed here. Because of the 
lower signal-to-noise ratio of the 2000 data not all the lines are 
found, so a comparison is only possible on the strongest lines. 
No lines appear to have statistically significant changes, possibly 
due to the poorer signal to noise in this data set. Examination 
of the spectra, how ever, reveals that several lines appear in the 
iBranduardi-Ravmont et al . ( 2001) spectrum but have fallen below a 
detectable level a year later. These lines occur at 0.510715±0.0003 
and 1.12631 ±0.001 keV and are identified as CaXIV and Fe XXIII 
respectively. 



6 DISCUSSION 

The analysis of the soft X-ray absorption lines of MCG-6-30-15 
have revealed some interesting properties of the warm absorber. 
Most importantly the column densities measured for O VII and 
O VIII lines are large enough, when considered with Fe I, to account 
for the drop in flux at ~0.7 keV. This means that the dominant ef- 
fect on the shape of the soft X-ray spectrum of MCG-6-30-15 is 
from a dusty warm absorber. Any additional complexity from rela- 
tivistically blurred emission lines is of low equivalent width. 

The column densities and turbulent velocities that were mea- 
sured for the other species agree with the predictions of the DWA 
model presented in Turner et aL, 12003.) . The measured equivalent 
widths of the identifi ed lines are larger tha n the predicted ones from 
the DWA models of iTumer et alj l l2003ft but only by a factor of 
~1.25 which may be accounted for by a slightly higher than So- 
lar abundance or higher than lOOkms"^ turbulent velocity in the 
warm absorbing gas than that assumed in the DWA models. 

A broad range of ionization states is found for several ele- 
ments, most notably oxygen and iron. This suggests the warm ab- 
sorber consists of man y zones with a range of ionization parameters 
as was assumed in lTumer et al . (2003) or possibly a conti nuum of 
ionization parameter as suggested bv ..Morales et alji200dh . 

The variability of various absorption lines suggests that the 
most lowly and highly ionized gas is varying as the continuum flux 
varies. The increase in equivalent widths of these lines with de- 
creasing flux is presumably from recombination of higher species 
in response to the drop in ionising flux. Using values for the re- 



combination timescale upper limits can be placed on the distance 
of these two zones from the source of ionizing flux. The low ion- 
ization zone gas is found to line within ~10 pc of the central source 
whereas the highly ionized zone must lie very close to the central 
engine at 10^^ cm, at a similar distance as the BLR. 
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